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I. INTRODUCTION 



In this paper we study the importance of nuclear finite size effects in inclusive muon capture reactions. The inclusive 
f*^ . muon capture process in nuclei 
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is very similar to neutrino scattering off nuclei 
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and is experimentally more accesible so it serves as a benchmark for testing theoretical models of the latter process. 

The motivation for this investigation comes from the results of a model (published in Q), which describes rather 
well inclusive 12 C(z^, /i~)A" and 12 C{v e ,e~)X cross sections near threshold and inclusive muon capture by nuclei. 
This approach, which is an extension of the quasi-elastic (QE) inclusive electron scattering model of Q , is based on a 
Local Fermi Gas (LFG), where the simplicity of the model makes it possible to include a great variety of effects into 
the reaction dynamics. The goal of this paper is to investigate whether finite nucleus effects that are not adressed by 
a LFG can affect significantly the results of the model in [l| . 

There already exist microscopic calculations of neutrino-nucleus reactions and muon capture that treat correctly 
the finite size of the system. However there are some dynamical issues that are implemented in a different fashion 
in a LFG model so a direct comparison of these models is not possible and we can not extract the effect of the 
t—j inclusion of nuclear structure details. For this reason, we perform a comparison of the LFG model with a extreme 
-h ' shell model (SM), that is single particle states in a Woods-Saxon (WS) potential, where the finite size effects can be 
easily recognized. 

In the second section we will introduce the model of p| as applied to inclusive nuclear muon capture and give 
some numerical results to be compared with experiment. Then we will introduce the extreme SM and a simplified 
LFG model. In the fourth section we will show the comparison between the two models. We will finish with some 
conclusions. 

More details on this issue can be found in 0. A more comprehensive analysis of the uncertainties in the model of 



can be found in 0. 



II. INCLUSIVE MUON CAPTURE IN NUCLEI 



The evaluation of the decay width for inclusive muon capture in finite nuclei proceeds in two steps. In the first one 
we evaluate the spin averaged decay width for a muon at rest in a Fermi sea of protons and neutrons T (p p , p n ) with 
N =/= Z. In the second step, we use the LFG approximation to go to finite nuclei and evaluate 

T = J d 3 r|0 ls (r)| 2 f(p p (r),p„(r)) (3) 

where <p± s (f) is the muon wave function in the Is state from where the capture takes place. It has been obtained by 
solving the Schrodinger equation with a Coulomb interaction taking into account the finite size of the nucleus and 
vacuum polarization. This approximation assumes a zero range of the interaction, that becomes highly accurate as 
long as the q dependence of the interaction is extremely weak for the /z— atom decay process. 

The spin averaged muon decay width is related to the imaginary part of the self-energy of a muon at rest in the 
medium. Further details and analytical expressions can be found in the Appendix of 0. 
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1939 


1310 


1311±8 


0.00 



Table I: Experimental and theoretical total muon capture widths for different nuclei. We quote two different theoretical results: i) 
Pauli+Q: obtained without including RPA correlations, but taking into account the value of Q; ii)RPA: the full calculation, including 
all nuclear effects. Experimental data (Exp) are a weighted average : r/cr 2 = E^/of, with i/ " 2 = E; V°i of thc results cited in 0. 
Finally, in the last column we show the relative discrepancies existing between the theoretical predictions given in the third column and 
the experimental data of the fourth column. (<5 rel r = (r Ex P - r Th ) /r Ex P) 

For kinematical reasons only the QE part of the W~ self-energy contributes to the muon decay. Thus, both the 
muon decay in the medium (Eq. ^) and the electroweak inclusive nuclear reactions (Eq. |2J) in the QE regime are 
sensitive to the same physical features. We can apply the same nuclear physics corrections to the above model as in 
neutrino scattering, that is Pauli blocking, RPA corrections and corrections to the energy balance, see yLj. The Is 
muon binding energy, Bjf > 0, is taken into account by replacing m M — > m M = m M — Bjf. 

In muon capture only very small nuclear excitation energies are explored, 0-25MeV, so the kinematical regime of 
the muon capture process is the worst possible for a LFG model of the nucleus. Nevertheless, the predictions of this 
model are in fairly good agreement with the experimental results as can be seen in Table [I] 

III. COMPARISON OF NON-CORRELATED MODELS 

In order to simplify the calculations in the comparison between the two models we make a static approximation 
and expand the single nucleon weak current in the nucleon momentum keeping terms up to order zero. 
In the SM we have to deal with an S matrix element of the kind 

S fi - -2mS(E f -Ei- ")^{f\ ^(-k')N) (4) 

where the weak current is modified to include the muon wave function 

J^(-k') = J d 3 re-* k '-%(r)0 ls (r). (5) 

Now the states |i) and (/| are nuclear states in a shell model. This states are single particle excitations of a N nucleon 
system in a WS potential and can be not only in the continuum, like in a Fermi Gas, but there can also be discrete 
excitations so we have now two contributions to the decay width. 

We now have to get the nuclear wave functions as solutions of the Schrodinger equation for a WS potential where 
the parameters of the potential are commonly fitted to the experimental energies of the valence shells or the charge 
radius. In the present case of muon capture we fit the experimental Q- value for the decay reaction. 

For thc LFG model we use simplified expression for the decay width, where we have not taken into account RPA 
correlations, but energy balance and Pauli blocking effects are implemented. The only inputs remaining to be fixed 
are the nuclear matter densities, that will be those provided by the wave functions of the WS potential. 

Up to now we have in both LFG and SM models the same physical features of Pauli blocking and correct energy 
balance, the only difference coming from the more refined treatment of the nuclear wave function in the SM case. In 
this way we can compare both SM and LFG models. 

IV. RESULTS 

In this section we present results for a set of closed-shell nuclei 12 C, 16 0, 40 Ca and 208 Pb. Fixing the experimental 
Q- value only makes one condition for fixing the several parameters of the WS potential, so wherever possible we set 
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Table II: Integrated width in units of 10 5 s~ for the different nuclei and Woods-Saxon potentials, compared with the LFG 
results using the corresponding charge densities. The discrete contribution of the shell model is shown in the first column. The 
column labeled % gives the relative diference between WS and LFG results. 
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Figure 1: The left figure shows proton and neutron densities of 1 C, for the several WS potentials used in this work. The 
right figure shows the differential SM width of 12 C to the continuum (left panels) compared to the LFG, and partial widths 
contributions to the discrete states (right panels) , as a function of the neutrino energy. 

the remaining parameters of the potential to values similar to those used in the literature. In our calculation we use 
different sets of parameters, denoted WS1, WS2 and WS3. In order to compare with the LFG, we use as input the 
proton and neutron densities obtained in the corresponding shell model. The values of the different parameters sets 
can be seen in Q . 

In Table HP we show results for the integrated inclusive widths for the four nuclei. We can sec that, in the case of 
WS1 and WS2, the LFG and WS results for 12 C are quite similar, differing only in ~ 2-3%. In the case of WS3 the 
differences are larger, around 14%. 

This can be understood in terms of the values of the parametrization WS3. For more attractive potentials the 
nucleus becomes more dense in the interior. For this reason, the WS3 neutron density turns out to be the smallest 
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Figure 2: The same as Fig. □ for 16 O. 

one, while the proton density is around 3/2 the neutron one. Therefore a proton near the Fermi surface can decay to 
a neutron above the neutron Fermi surface with an energy decrement. This is an unrealistic situation, since precisely 
in this case the neutrons are less bound than protons in the SM, and therefore lie at higher energies. Hence the LFG 
results are worse for very different neutron and proton densities. Another argument to disregard this case is the well 
known property of closed (sub-)shell nuclei such as 12 C, for which the neutron and proton densities should be similar. 

In Fig. ^ we compare the SM results for the differential width to the continuum with the LFG distribution for the 
different WS parameters (left panels). The shapes of both distributions are completely different. The partial widths 
to the discrete states are shown in the right panels of Fig. 2] Considering these differences in shape between the LFG 
and the SM, it is a very notable result that the integrated widths (adding the discrete states) take similar values in 
both models as was shown in Table ITT1 

In the case of 16 the integrated widths computed in the LFG are also very close, ~ 4-6%, to the SM results 
with the potentials WS1 and WS3 (see Table ITTI) . The worse results are obtained for the WS2 parameterization; the 
corresponding width is 24% of the SM one. This can also been understood in terms of what was said for the case of 
12 C above, by looking at the 16 densities shown in Fig. [3 

The LFG results improve when the mass of the nucleus increases as in the case of the nucleus Ca. In fact, from 
Table ITTI we see that for this nucleus the LFG integrated width is lower than 8% for all cases. This improvement was 
expected because the Fermi gas description of the nucleus should work better for heavier nuclei. 

For the closed-shell heavy nucleus 208 Pb we present in Table ITT) integrated widths only for two sets of potential 
parameters, WS1 and WS2. In both cases the LFG results are close to the SM ones. 

V. CONCLUSIONS 

In this paper we have estimated the magnitude of the finite nucleus effects on inclusive muon capture, aiming at 
quantifying the uncertainty of the LFG results of jl|. 

We have focused on a simple shell model without nuclear correlations, but that contains the relevant information 
about the finite nuclear structure, and we have compared it with the uncorrelated LFG using the same input. As 
expected, the neutrino spectrum is very different in the two models, in particular the LFG cannot account for the 
resonances and discrete states. However, in the case of the lighter nuclei, 12 C and 16 0, the SM and LFG results for 
the integrated width are close — within 3-6% — for WS parameters with similar neutron and proton densities. For 
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Figure 4: The same as Fig.^Jfor 208 Pb. 



the medium and heavy nuclei, 40 Ca and 208 Pb, the integrated widths are always very close, within 1-7%. The final 
neutrino spectra of the LFG become more similar to the SM, including the discrete part, for heavier nuclei. 
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